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Abstract 



We first obtain the most general and compact parametrization of the unitary transformation 
diagonalizing any 3x3 hermitian matrix H, as a function of its elements and eigenvalues. We 
then study a special class of fermion mass matrices, defined by the requirement that all of the 
diagonalizing unitary matrices (in the up, down, charged lepton and neutrino sectors) contain 
at least one mixing angle much smaller than the other two. Our new parametrization allows us 
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^ to quickly extract information on the patterns and predictions emerging from this scheme. In 

particular we find that the phase difference between two elements of the two mass matrices (of the 
sector in question) controls the generic size of one of the observable fermion mixing angles: i.e. just 
fixing that particular phase difference will "predict" the generic value of one of the mixing angles, 
irrespective of the value of anything else. 
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I. INTRODUCTION 



In the absence of flavor symmetries, the Yukawa couplings between the Standard Model 
(SM) fermions and the Higgs field are in general complex arbitrary matrices which, after 
Electroweak Symmetry Breaking (EWSB), become the mass matrices of the quarks and 
charged leptons. In the case of neutrinos, the mass matrix will be in general complex 
symmetric. All these matrices contain more parameters than physical observables and an 
explicit computation of these observables (fermion masses and mixings) in terms of the 
original matrix elements can be quite cumbersome in general. Indeed this would require us 
to solve a 3 x 3 eigenvalue problem for each fermion matrix, and then compose the unitary 
transformations (formed with the calculated eigenvectors) of the Up and Down quark sectors 
and then also of the charged lepton and neutrino sectors 1 . We observe however that it 
may be useful to address the question not as an eigenvalue problem, but as an eigenvector 
problem, treating the eigenvalues as input parameters and not as output. The first reason 
for this is that except for the neutrino sector, all the mass eigenvalues are quite well known. 
But the main point we make is that by keeping explicitly the mass eigenvalues as input 
parameters, the eigenvector solutions of each mass matrix become surprisingly simple and 
can be written as compact functions of both the mass matrix elements and the eigenvalues 
(the fermion masses). Such a parametrization of the mixing matrices, directly in terms of 
the original mass parameters and the fermion masses might prove to be useful in the studies 
aiming to explain the observed flavor structure of the SM by way of symmetries or textures 
or patterns [IHI] at the level of the fermion mass matrices. 

It is true, though, that to proceed we need to work with hermitian matrices, but it is 
always possible to render the quark mass matrices hermitian in the Standard Model without 
loss of generality [51 E] • The procedure to obtain hermitian matrices is quite standard, and 
it involves either working with the hermitian matrix Hi = GG\ where G is the original 
fermion mass matrix, or using its polar decomposition, i.e. solving G = H 2 Q G , where Q a 
is a unitary matrix converting the general complex matrix G into a positive semi-definite 
hermitian matrix H 2 (if G is invertible, with distinct non-vanishing eigenvalues, then H 2 is 
positive definite and therefore Q G is unique). 

We will first present our parametrization (in fact 9 different types) for the unitary matrix 
W which diagonalizes a general hermitian matrix H in the most compact way possible. Then, 

1 Assuming three families of neutrinos, although additional neutrino flavors are possible (sterile neutrinos). 
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to start taking advantage of it, we then propose a simple and mildly constraining ansatz for 
the flavor structure of the SM fermion sector. It assumes that the unitary transformations 
diagonalizing H u and Hd can each be decomposed as only two rotations, instead of three. 
The idea is to assume that the third rotation angle is zero (or much smaller than the other 
two) and therefore one of the entries of each transformation matrix W u and Wd will be zero 
or close to zero. We call this setup the two-angle ansatz and we will concentrate in only 
one of the many possible cases. Our parametrization allows us to quickly obtain very simple 
dependences of the fermion mixing matrices Vqkm arid Vpmns [Zj m terms of the masses 
and the original mass matrix elements. We can thus study easily the interesting properties 
of this type of ansatz as well as the consequences it has in the original mass matrices, in both 
the quark and the lepton sectors. A particularly interesting observation is that the specific 
value of some elements of H u and Hd has no effect (or very mild effect) on the observed 
values of masses and mixings. 



II. MIXING MATRIX PARAMETRIZATION 



As explained before, we are going to concentrate on hermitian matrices with the as- 
sumption that the fermion mass matrices are either hermitian or that one can construct a 
hermitian matrix out of them. All the results presented in this work are valid for general 
hermitian matrices, but for simplicity we will only consider the case of positive definite 
hermitian matrices. Let H be one such matrix: 

(1 



H 



x 



X 

a 
If 



g\ 
b 

a ) 



(1) 



with eigenvalues Ai, A2 and A3. A compact parametrization of the unitary matrix W which 
diagonalizes it, is 2 

gb* — x(a — X 2 ) 

1 2 



/ (q-A 1 )(a-A 1 )-|b| s 



W 



g*b — x*(a — Ai) 



Ax 

x*b* -g*\a- Ai) 
Ai 



N 2 

(7 - AgXg - A 2 ) 



xb — g(a — A 3 



\ 



\g\ 



N 2 

xg* - 6* (7 - A 2 ) 
N 2 



A 3 

x*g - 6(7 - A 3 ) 



, ^3 

(7 - A 3 )(a - A 3 ) - \x\ 
A 3 



(2) 



2 There are 9 different ways of parametrizing it, depending on the choice of which diagonal mass parameter 
(7, a or a ) is explicitly absent in each vector column (see Appendix A for details). In the parametrization 
shown here, the rotation matrix W has the correct limit when the off-diagonal entries in the original mass 
matrices H are set to zero, avoiding (apparent) divergences in this limit. 
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After some algebra the normalization parameters are found to have the simple form 

Nl = (A 3 - Ai)(A 2 - Ai) [(a - Ai)(a - Ai) - |6| 2 ] , (3) 
N 2 2 = (A 3 - A 2 )(A 2 - AO [{a - A 2 )(A 2 - 7) + \g\ 2 ] , (4) 
Nl = (A 3 - A 2 )(A 3 - Ai) [(A 3 - 7) (A 3 - a) - M 2 ] (5) 

The surprisingly simple and compact form of this parametrization might make it suitable to 
treat flavor models keeping always an explicit dependence on all the matrix elements of the 
hermitian mass matrices. Of course, if the three eigenvalues Ai, A 2 and A 3 are fixed, there 
must be three constraint equations on the elements of the matrix H. These equations are 
found from the three invariants Tr(H), Tr(H 2 ) and Det(H): 

Tr{H) = a + a + 7 = Ai + A 2 + A 3 , (6) 
Tr(H 2 ) = 2(\x\ 2 + \b\ 2 + \g\ 2 ) + a 2 + a 2 + 7 2 = A 2 + A 2 + A 2 (7) 
Det(H) = 7 (cra - \b\ 2 ) - a\x\ 2 - a\g\ 2 + 2Re(bxg*) = AiA 2 A 3 (8) 

By choosing g as an independent variable, it is possible to rewrite these constraint relations 
on the rest of variables as 3 

a = Ai + A 2 + A 3 - a - 7, (9) 
2 = (7 ~ AQ(A 2 - 7) (A 3 - 7) ~ \g\\a ~ 7) + 2Re(bxg*) 
X (a - 7) 

2 _ (a - Ai)(a - A 2 )(A 3 - a) + \g\ 2 (a - a) - 2Re(bxg*) 

(a -7) 

The interesting thing of this notation is that the constraint formulae on x and b actually 
become algebraic solutions for both x and b when the term Re(bxg*) vanishes identically. 
In particular, this is the case when one deals with mass matrices with texture zeroes in the 
off-diagonal elements. 



III. FLAVOR IN THE TWO- ANGLE ANSATZ 



Equipped with an exact and simple parametrization of the fermion mixing matrix in 
both Up and Down sectors (or charged lepton and neutrino sectors), we look for economical 
patterns among the mixing matrices by following a bottom-up approach in the hope that it 
might be complementary to more top-down approaches such as imposing flavor symmetries 

3 The same type of relations can be written when choosing x or b as the independent variable. 
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or texture zeroes mass matrices (see [IHI] as well as the probably incomplete surveys of 
[HI E])- One avenue is to find a similar ansatz for the flavor structure of the mixing matrices 
in both quark and leptonic sectors. Such possibility exists in the sense that in both sectors, 
the mixing elements (Vckm) 13 and (Vpmns) 13 are known to be small. This feature is very 
interesting and it is known to lead to a simple parametrization of the mixing matrix. Note 
that in the limit V13 = we have two additional conditions, since V13 is a complex number. 
This means that in this limit, the mixing matrix will have only two independent parameters 
instead of four. We can choose these four independent parameters to be | V12 1 , | V23 1 , | V^3 1 
and I V21 1 [ID]- In the limit of V13 = 0, we have the extra constraint | V21 1 = |Vi2|\/l — |V2 3 | 2 . 
This means that the whole mixing matrix is described by | V12 1 and | V23 1 . Note also that in 
this limit, there is no CP violation a la Dirac in both quark and leptonic sectors. 

In what follows, the subscript stands for the values of the mixing matrix elements in 
the limit Vi 3 = 0. Since V13 is known to be very small, we believe that the zeroes values 
of the mixing matrix elements are not far from their measured values. For the quarks, we 
have: 
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Note that this zero-order Vq KM can be decomposed as a product of two rotations, namely 
one is purely the Cabbibo angle and the other one is purely made out of beauty namely 
\V cb \: 

VSkm = VbV c (13) 
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In the leptonic sector we have: 
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/ 



which can be decomposed also as a product of two rotations, namely one is purely solar and 
the other one is purely atmospheric: 
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(16) 
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The diagonal phase matrix P which contains the Majorana phases is defined as: 

/l \ 



P = 











V 



(17) 
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This structure for the physical fermion mixing matrices, decomposed mainly into just two 
rotations is quite suggestive and we will use this observation as the starting point of our anal- 
ysis. In this context we ask ourselves how many economical possibilities one has to restore 
minimally and fully the mixing elements and CP violation. We start with the structure 



(19) 



and then to be as general as possible while still keeping the original motivation we consider 
all mixing patterns emerging from the original structure, but with only one zero in each 
mixing matrix. Now, we establish all the corrected mixing patterns for both Up and Down 
quark sectors: 
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where X stands for a non-zero value and Cor for a corrected originally zero mixing 
matrix element. These mixing matrices with one texture zero in them can be decomposed 
themselves into only two rotations instead of three. We call this the two-angle ansatz, in 
which both the Up and Down quark sectors are diagonalized by unitary transformations 
containing just two angles (i.e. having one vanishing element). There are obviously many 
possibilities for this ansatz but in particular there are 16 cases such that one can recover in 
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a specific limit the case Vj.3 = 0. In this work we are going to focus only on one specific 
example of this type of ansatz, although a full case by case study is underway. We feel that 
the main features of this ansatz do reveal themselves in the example studied here, and we 
prefer to continue elsewhere a more systematic exploration. 



Notation 

Before we proceed further, we will set up our notation for all the matrices in both up and 
down quark sectors as well as charged lepton and neutrino sectors. We will define the 
Hermitian matrices H u , Hd, Hi and H v defined by: 



(22) 
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with eigenvalues m u , m c , m t 
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(23) 



with eigenvalues m e , m^, m T , Ai, A 2 and A 3 . 

In the text, we will also need to refer to the phases of the off-diagonal terms, denoted as 



arg(# M )i2 = 5 X 
arg(if d )i2 = 6. 



13 - S g , 



arg (H u 
arg(F d )i 3 = 5 h 



arg (H u ) 23 = S b , 
arg(# d ) 23 = 5 f , 



(24) 
(25) 



and similarly for the leptons. The matrices diagonalizing these mass matrices will be denoted 
respectively as W u , Wd, Wi and W u . 



A. THE CASE (13-13): i.e. (W u ) 13 = (W d ) 13 = and/or (Wi) 13 = (W,) 13 = 



We will now consider the mass matrices from Eqs. (22 ) and (23 ) with the extra constraints 
of (W u )i3 = and (Wd)i3 = in the quark sector and (Wj)i 3 = and (W„)i3 = in the 
lepton sector. 4 



4 In the quark sector, this limit was shown to lead to acceptable patterns in the limit p = and 7 = [TT] 
and similar limits were considered in the lepton sector in |12) 
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1. The quark sector 



For example in the down quark sector imposing (Wd)i3 = corresponds to the require- 
ment yf — h(/3 — m b ) = . After a short computation, we obtain simpler relations among 
the elements of the mass matrix: 
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\h\ 2 
l/l 2 



(p - m d )(m s - p){m b - (3) 
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(m b - 0)(m b - d) 



S v - S h + 5 



f 



71 



(26) 



These relations will simplify significantly the general form of the Down quark mixing matrix 
Wd as well as the up quark mixing matrix W u (see Eq. (|2| and Appendix B for details). 
It is then straightforward to show that the quark mixing matrix takes the form 



v = wlw, = Pi v r 



p 

U " CKM d 



(27) 



where P u and Pd are unphysical diagonal phase matrices given below and 
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The complex rotation parameters C, S, c 7 and c p are given by 
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(i.e. \C\= \V tb \) (29) 

(30) 

(31) 
(32) 



Note that S; = a/1 - c? and that \C\ 2 + \S\ 2 = 1. The CP phase 6 and the two 
unphysical phase matrices P u and Pd are given respectively by 6 = 5 X — 5 y — 5c , 
P u = diag(l, e i5 % e i( ^ +5s) ) and P d = diag{\, e i& \ e i(5a;+<5s) ), with 5 { = arg (z). 
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The previous form of the mixing matrix implies the following exact relations for the quark 
sector 



\Vub\ h - m u \V td \ [p - m d 



\V cb \ V m c-7' \Vu\ \m s -p 



(33) 



and 



\V ub \ (l - m u )(m s - m d ) \V cb \ \{m c - i)(m a - m d ) 



\V t d\ V ( mc _ m u)(p ~ m d ) ' \V ts \ V (m c - m u )(m s - p) 



(34) 



which really correspond to two independent constraints (i.e. the third and fourth relations 
can be obtained using the first two along with the unitary constraints). 

As can be seen only two original mass matrix elements 7 (from H u ) and p (from H d ) 
appear explicitly in these last four relations showing the first effect of the (13 — 13) ansatz, 
i.e. linking each of the previous ratios to a different quark mass matrix, and in particular 
to the first diagonal elements of each mass matrix. Once we fix these two elements to fit 
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the experimental value of the ratios given in Eqs. (33) and (34), we will still have 6 free 
parameters (including 4 phases) to fit the rest of the data (i.e we need to fit two more scalar 
observables from experiment, which for example can be taken to be the absolute value of 
\Vtb\ and the value of the CP phase of the fermion mixing matrix). We can choose to use a, d 
and the phases of x, g, f and b as the free parameters. It is also useful to define the phase 
combinations: 

(35) 
(36) 
(37) 

which have the constraint 

A 23 +A 12 = A 13 . (38) 

due to the (13 — 13) ansatz imposition. 

Of course, we have more than enough free parameters to fit the two remaining observables, 
but as we will shortly see, of the 4 free phases, only two phase differences can be relevant, 
and the 2 real parameters a and d turn out to be statistically irrelevant if the two phases are 
properly chosen. In other words, in most of the parameter space of the six free parameters 
needed to produce the two remaining physical observables, 4 directions are more or less 
irrelevant, with two phase differences being the two parameters required to obtain a good 
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FIG. 1: Distribution of | Vtb{a, d, &23)\ with random a, d and A23 with respect to A23 (left panel), 
where a and d are two diagonal elements of the quark mass matrices H u and Hj, and A23 is the 
difference between the phases of the ^(23) elements of these matrices. On the right panel, we present 
the distribution of the Jarlskog invariant J(a, d, A12) with random a,d and A12, for A23 = 0. On 
the two panels, the black curves represent the 25% and 75% quantiles of the distribution for fixed 
A23 (left) and A12 (right). In other words 50% of the random points lie between the curves, with 
25% above it and 25% below it. 



experimental fit. This situation is somewhat surprising because the amount of "useful" free 
parameters is less than the total number of free parameters. Let's see how this is played 
out. 



Eq. (33) shows that 7 = 0{m u ) and p = 0{m ( i) are required in order to obtain a good 
fit with experimental data [13J. This has interesting implications for V t b since it means that 
a + a = m t + m c + 0{m u ) and d + (3 = nib + m s + 0(md) after using the trace identity of 
H u and Hd- We can therefore write 



V t b(a,d,A 



23 I 



y/(m t - a)(m b - d) e lA23 + y/(a- m c )(d - m s ) 



(39) 



y/(m t - m c ){m b - m s ) 

It turns out that when the phase A23 is small (modulo 2n), statistically we find that 
\Vtb{a,d,0)\ ~ 1 for any randomly chosen value of a and d. In fact, the generic value of 
\V t b(a, d, A 2 3)| is very much correlated with the value of A23, with little dependence on the 
values of the other two variables, at least for small enough A23. This is shown in figure [TJ 
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Regions where | V t b I > 0.995 for various A 2 3 




50 100 150 



a (GeV) 

FIG. 2: Regions in the plane (a,d) where \Vtb\ > 0.995, for three different values of A23. The 
dashed curve is a contour of the experimental fit for Vtb — 0.9991 in the case A23 = 0. One sees 
quite clearly that as soon as A23 is increased, the parameter space favorable for large values of Vtb 
shrinks dramatically. 

where we plot the distribution of | V t b(a, d, A 2 3)| with respect to A23 for randomly chosen 
values of a, d and A23. It is apparent that there is a clear correlation between the value 
of the phase A23 and the value of \V t b\- The two black curves correspond to the 25% and 
75% quantiles of the distribution of \V t b\ for a given value of A23 (i.e. 50% of the randomly 
generated points lie between the two curves, with 25% above them and 25% below). 

All this shows that the generic value of V t b in this ansatz is actually governed by the 
specific value of the phase A23, with very mild dependence on the other 2 parameters a and 
d. It is not clear though, that demanding A 2 3 ~ will be enough to fit the observed data in 
the quark sector since V/ b lt ~ 0.999, i.e it is quite close to 1. The random scan of Figure [l] 
does not show the distribution of V t b for A23 = 0, but once we fix A23 we just have two free 
parameters, a and d. In Figure [2] we show the complete allowed phase space for a and d, 
which are subject to the experimental constraints m u < a < m t and rrid < d < rrib . The 
different shaded regions are the points where Vtb > 0.995, for three different values of the 
phase difference A23. When this last one is zero, the region is quite large, and it decreases 
very fast as the phase is increased. Because the area of parameter space can be quite large 
we may still say that the value of \Vtb\ consistent with experimental data could be considered 
"generic" as long as the phase difference A 2 3 is vanishingly small. 

Finally, one can also compute quite easily the Jarlskog invariant [H] in this context. 
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From Eq. (28), and using for example the definition J = Im(V cb V t * b V* d Vtd) , it is easy to see 



that it will have the form: 



J = CpCy SpS-y \C\ \S\ 2 sin 9 



(40) 



where q, C and S were given in Eqs. (29), (30), (31) and (32), and where 9 = S x — 5 y — 5c 
with 5c = argC. 

Before analyzing in more detail the dependence on the original mass matrix elements, it 
is interesting to relate the phase 9 not just to the Jarlskog invariant but also to the angles 
of the unitarity triangle <pi — (3 , 2 = a and 3 = 7. These angles are defined as: 
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and the relation between them and the Jarlskog J is: 



(41) 
(42) 
(43) 
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In the case of our ansatz, we can rewrite our J as 
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For example, from the first identities of Eq. (44) and (45) one sees that we must have 
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c p c 7 sin 6* 



IK 



ud 



(46) 



Since c p c 7 ~ \Vud\, h follows from the above equation that a ~ 9. The experimental 
constraints on a are such that a = (88^5)° [13], which basically means that the phase 9 is 
constrainted to be 9 ~ 7r/2. 

A more revealing way to see what this means for the original elements of the mass 
matrices, we can actually rewrite the Jarlskog in a more explicit way as 



J = \V ub \\Vtd\cpC^ 



- \V ub \\Vtd\cpC~f 
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(47) 
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where the approximation of the second line comes from assuming that a + a = m t + m c + 
0{m u ) and d + (3 = m b + m s + 0(m,i) . 

We showed earlier that the imposition of A 2 3 = (Sf — 5 b ) ~ gives a nice statistical 
reason for the large value of \V t b\- With that extra condition, one actually has A 12 ~ 6, 
which then means that A 12 ~ it/ 2 due to experimental bounds. The experimental best fit 
value of the Jarlskog invariant is \Jfu\ = 3.05 x 10~ 5 |I3] ; in the right panel of Figure [l] we 
present a scan of values of the function of J shown in Eq. (48), for random values of a, d 
and A12 (assuming A23 = 0). It is apparent that the observed value of J can be obtained 
quite generically when A12 = tv/2, in a quite insensitive way to the specific values of a and 
d. 

It is quite suggestive that some specific phase differences between the elements (-£/«) 23 = 
/ and (Hd)23 = b , and between the elements (H u )i 2 = y and (Hdjn = x and between 
the elements (H u )i 3 = h and (H d )i 3 = g , given by: 

A 23 = S f -5 b = (49) 
Aia = 5 y -8 x = I (50) 

A13 = S h - S 9 = I (51) 

do lead to good generic values of both V t b and J. We are left with four parameters a, d and 
two combinations of phases independent of A12 and A23, all from the original mass matrices, 
which do not seem to play any important role in obtaining a good overall fit in the quark 
sector. 

Before we finish this subsection on the quark sector, we would like to point out that since 
a and d can take almost any value (inside their allowed range), one might actually get very 
close to a symmetric limit, namely a (2 <H- 3) family symmetry, relating second and third 
families (see [9] for implementations mostly in the lepton sector). 

Forgetting for a moment our ansatz but assuming the (2 -H- 3) limit for both up and 
down quark mass matrices, one is then guaranteed to have vanishing elements (Wu)i3 and 
{Wd)i3 (i.e. we recover our ansatz). Moreover some of the elements of the mass matrix are 
subject to the constraints 

x = g y = h (52) 

a = a (3 = d (53) 

6 b = TT = 6 f . (54) 



From these equations, we must have A 2 3 = 5/ — 5 b = and using Eq. (29) it is easy to see 
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that \Vtb\ = 1 exactly, which means also that \V u b\ = \V c b\ = \Vtd\ = \V ts \ = 0. Since the 
fitted values of \V u b\, |V^|, \V t d\ and \V ts \ are at most 10~ 2 and much smaller than V t b, then 
it seems plausible that a small deviation from this symmetric limit can easily restore the 
experimental values of these Vy's. 

To obtain a correct experimental fit, one would also have to require the phases of x and y 
to be separated by tt/2 (as remarked earlier), and moreover the values of p and 7 will have 
to be chosen so as to obtain the correct ratios and 1^4 (see Eq. (33)). 

Although we imposed a flavor ansatz mainly for empirical reasons and simplicity, it is 
interesting that the patterns emerging from it do actually lead to a possible type of symmetry. 
Of course we just considered one of the possible cases of our scheme, and others cases might 
lead to different patterns or reveal some other feature and for this a more systematic study 
is required, some of it being already underway. 



2. The lepton sector 



In the lepton sector, the mixing matrix takes the same form as for the quarks in the 
(13 — 13) ansatz, i.e. 



r I ' ' V ± I y PMNS 

where Pi is an unphysical diagonal phase matrix given below and 

/ CyC p i + S p >Sy\C'\e ie ' — CyS p ' + SyC p '\C'\e t9 ' Sy\S'\ \ 



(55) 



-SjiCpi + Cy Spf I C I e %9 ' SyS p i + Cy c p > \ C | e 10 ' C^i I S' 



V 



-s p >\S'\ 



-Cp>\S'\ 



T 

\C'\e- id ' J 

where P v is the Majorana phase matrix given below and the complex rotation parameters 
C',S', cy and c p i are given by 



(56) 



C 



S' = e ld » 



v /(m r -a')(A 3 -rf / ) eW-W + = «')(A 3 - P')) 

v/(2m T -a'- a')(2A 3 - (3' - d') 
xi U{m T - a')(A 3 - W) - eWr-Vy/fc - a')(A 3 - d'j) 

id, \ J 



v/(2m T -a' - a')(2\ 3 - J' - d') 



Cy 



m p 


-i 


m p 


- m e 


'A 2 - 


-P' 



A2 — Ai 



(i.e. \C'\ = \V T3 \) (57) 

(58) 
(59) 
(60) 
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Note that sj = 1 — cf and that \C'\ 2 + |S"| 2 = 1. The CP phase 9' and the Majorana phase 
matrix P u , as well as the unphysical phase matrix P; are given respectively by 9 = 5 x r — 5 y > — 
S c , , P u = diag(l, e iS y\e i{ - 5 -' +6 c)) and Pi = diag(l, e iS *' , e i( V +<5 c)) ; wnere §. = a rg(i). 
Again, the previous form of the mixing matrix implies the relations: 



\v e3 




lj'-m e 








m M - i ' 


^r2 



y-Ai 



(61) 



and 



|Ves| _ /(7 , -^e)(A 2 -A 1 ) |V^ 3 | _ /(m M -7')(A 2 -Ai) 



|Ki| V ("V _ m e)(p' - Ai)' \V r2 \ V (w„-m e )(A 2 -//)' 



(62) 



with again only two parameters from the original mass matrices separately controlling each 
ratio. 

But, the lepton case is different, obviously, because the lepton mixing matrix contains 
two large angles and also because the neutrino mass structure is unknown and might not 
be as hierarchical as in the quark sector (see for example [15] for the latest global fits 
coming from neutrino oscillation experiments). In the (13-13) ansatz we see that we must 
still enforce the charged lepton mass matrix to have a very small first diagonal element, i.e 
7' = 0(m e ) in order to obtain a small value for the ratio | 1 /| V^3 1 . This will imply that 
a' + a' = + m T + 0(m e ) (from the trace identity of Hi) which will simplify the functional 



form of V T z of Eq. (57). Before doing so, we need to look also into the required value of p' to 
obtain a correct fit for 1^4 ~ 1/ y/2 in the approximate tri-bimaximal scheme (TBM) [TB] . 
This leads to p' ~ A 2 + 2Ai, and we can now use this relation to obtain the simplified form 
of V r3 : 

(^(mr - a')(D 13 + \! - d')e iA '™ + y/(a' - m„)(Ai + d')) 

V T3 (a', d', AL, Ai) = ^ , '- (63) 

v /(m T -m M )( J D 13 + 2A 1 ) 

where D\ 3 = (A3 — Ai) is fixed by the measured atmospheric neutrino mass difference 5 . 

We now have 7 free parameters (including 4 phases), since one of the neutrino masses is 
unknown, and we choose it to be Ai = m 2 . As seen in Figure [3j when mi < 0.1 eV, the 
correlation between the phase A' 23 and the value of \V r3 \ is basically the same as it was in 
the quark sector, and this is basically due to the fact that this feature happens for masses 

5 Although we have been writing Ai, A2 and A3 for the neutrino eigenvalues, these could in fact correspond 
to the squared physical masses if we consider the neutrino mass matrix squared as our starting point. In 
this case, D 13 = (Am 2 ) atm . 
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Random Distribution of | V t3 (a ', d', A '23, m{) | Random Distribution of J (a', d\ A '12, »«i)/ I \ 




FIG. 3: Distribution of \V T 3(a' ,d' , A' 23 ,mi)\ with random a', d', A 23 and mi (denned as the 
mass of the lightest neutrino) with respect to A 23 (left panel). On the right panel, we present 
the distribution of the Jarlskog invariant J(a', d', A^ 2 , mi) normalized to the value of \V e s\, with 
random a',d', A' 12 and mi, for fixed A 23 = tt/2. On the two panels, the black curves represent the 
25% and 75% quantiles of the distribution for fixed A' 23 (left) and A' 12 (right). In other words 50% 
of the random points lie between the curves, with 25% above it and 25% below it. 

which are hierarchical (and in fact, for larger values of m\ the correlation between the phase 
A' 23 and the value of | 1 starts to wash out). It is quite suggestive, that if m\ < 0.1 eV and 
A 23 = 7r/2, then we obtain a generic size of jV^] ~ l/y/2, as shown in Figure [3} Because 
7' is small, cy ~ 1 and we will also have |V^ 3 | ~ l/\/2. These two values of |V^ 3 | and \V T s\ 
are consistent with the tri-bimaximal scheme (TBM) where | V^^^^ | = IK^" 8 ^ = l/v2 and 
so we find that in the lepton case the preferred value of the phase combination A' 23 is |, 
which statistically predicts a large value for both the mixing angles V T % and V^ 3 . It is also 
interesting to note that the charged fermion mass matrix must have the same structure as 
the up and down quark mass matrices, namely that the first diagonal element is of the order 
of the lightest eigenvalue. 

Since the mixing element V T \ is not small, we also conclude that the neutrino mass 
structure must be different from the other three matrices. This could be due to the fact 
that the eigenvalues are not as hierarchical as for the charged fermions. As in the quark 
case, it is interesting to note that further imposing some flavor symmetry might simplify the 
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relations obtained in this ansatz. In particular one could ask what is the effect of imposing 
the (2 ■<->• 3) family symmetry (see for example [H]) only in the neutrino sector. In our case 
this is a natural question to ask since this symmetry imposes automatically the element 
(W^)i3 to be zero, which is the defining condition of our ansatz. But it also imposes the 
following constraints on the mass matrix elements 



arg (H v ) 23 = 7T, 



i.e. 
i.e. 
i.e. 



V = h' 
(3' = d' 

6f> = 7T 



(64) 
(65) 
(66) 



This will have interesting consequences on the lepton mixing matrix because in this case we 
will have 

J - y/(m T - a') e- iS »' 



t3 



a 



i 



m, 



V2 



y/(m T 



+ C(m e /m M ) 



(67) 



where we have assumed that 7' = 0(m e ). Now one can see that when 5y = vr/2, independent 
of anything else, we will have \V t3 \ = 1/V2, as well as \S'\ = l/v2. Thus the 2^3 
symmetry, along with the requirement of the phase of b' to be tt/2, forces the appearance 
of large mixing angle in the Vpmns matrix, and we remind the reader that we are not 
considering a diagonal charged lepton mass matrix, although all our matrices are taken in 
the (13 — 13) ansatz. The form of the Vpmns matrix under these assumptions is 

c Y /V2 P v (68) 
-c pl /V2 e- i6 '/V2j 

and where 6' is the Dirac CP phase 



/ CyCpi + S p /Sye l8 ' /a/2 —CyS p i + SyC p / e l6 / \f2 



PMNS 



sySpi + cyCpie 16 ' /\/2 



where as before, <v 



and Cpi 



"V y m^—vrif, P y A2— Ai 

given by 9' = 5 X > — S y > — 5c and the Majorana phase matrix P v is given by 



/I 






e i5 v' 








\0 



(69) 



The phase 8c> is given by sin 5, 



c> 



e i(6 x ,+8 c ,) J 

f 0{m e /mp). 



Finally, in the more general case where one does not consider (2 <f->- 3) symmetry in the 
neutrino sector, the leptonic Jarlskog invariant J is computed to be: 



J = |^3||Ki||K3|c;<sin0' 



|K3l|Ki|c'/ 7 - 



^(m T -a')(rn3-d') sin (A^ 3 +A' 12 ) + y/ (m T -a')(m 3 -(3') sin A' 1; 
y/(2m T -a' -a')(2m 3 -P' -d') 
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(70) 
(71) 



For the preferred value of A' 23 = it/2, such that statistically it is more favorable to obtain a 
V T 3 mixing angle close to 1/ \/2, we obtain 

("\/ '(m T -a')(m 3 -d') COS A' 12 + y/(m T -a')(m 3 -P') Sm A[ 
y/ (2m T -a' -a')(2m 3 -p" -d') 



J = \V e3 \\V Tl \c> p c^ ______ L (72) 



This suggests that now, the value A' 12 ~ 7r/4 will maximize the possible value of J (normal- 
ized to | V e s | ) , and this is confirmed in the random scan shown in Figure 13} 



IV. DISCUSSION AND OUTLOOK 



In this work, we started by obtaining the simplest parametrization for the diagonalization 
of a 3 x 3 hermitian matrix, as a function of both its matrix elements and its eigenvalues. 
Since the masses of fermions in the SM are well known, the problem to attack is not an 
eigenvalue problem, but an eigenvector problem. In other words, we are able to obtain the 
eigenvectors of a hermitian matrix in an algebraically compact form because we treated the 
eigenvalues as known parameters, instead of unknown. With this parametrization in hand, 
we wanted then to show how it could simplify and make quite transparent the analysis and 
study of some flavor schemes. To this end we defined a new flavor scheme which imposes some 
(arguably obscure) constraints on the flavor structure in the fermion sector. The constraint 
imposed is the requirement of the vanishing of one of the mixing angles of the diagonalization 
matrix of each hermitian fermion mass matrix of the SM. We called this the "two-angle" 
ansatz (i.e. out of three angles to diagonalize a hermitian matrix, we consider the family 
of hermitian matrices diagonalized by only two angles), and it is bottom-up motivated, i.e. 
it is inspired on the observed structure of both Vckm and Vpmns matrices. Nevertheless 
we also observe that in flavor models where one requires a full symmetry between two of 
the three families, like the y, — t symmetry models, one of the consequences is precisely 
the vanishing of one of the mixing angles. 6 But we mainly decided to keep the bottom- 
up approach motivation and study the patterns emerging from our "two-angle" ansatz and 
focused on just one possible implementation, the (13 — 13) case, in which the (W / ) 13 entries of 
all the diagonalizing matrices happen to vanish. Using our parametrization, one can actually 
write both Vckm and Vpmns in terms of the original matrix elements and eigenvalues. In 
particular we found a peculiar dependence of one of the mixing angles with the model 

6 In the case of a complex symmetric neutrino mass matrix, the authors of [IT studied the symmetry 
conditions required for the vanishing of one of the mixing angles of the rotation matrix 
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parameters. We observed that by fixing the phase difference A23 between two parameters 
of the up and down mass matrices (or the neutrino and charged lepton mass matrices), 
the mixing angle V t b (or V T z) can be "predicted" in a statistical sense, i.e. if one makes a 
random scan allowing the remaining free parameters to take any possible value, including 
the eigenvalues, one obtains a narrow distribution for the value of \V t b\ (or |V^ 3 |), and the 
central value is a monotonic function of the phase difference A23. This can be seen by 



rewriting the general formula for V t b in the (13 — 13) ansatz (see Eq. (29)) as 



\V t b\ = F(a, a, d, (3, A 23 ) = Abs 



( ( yj(l - a)(l - d) + y/(l-a)(l-P)\ ^ 

y / (2-a-5)(2-/3- J) 



(73) 



where a = a/m t , a = a/m t , d = d/rrib and /3 = (3/mb, and a, a and d, /3 are diagonal 
elements of the up and down mass matrices respectively. If the mass eigenvalues m, = A"' d 
are unconstrained, the only constraints on the parameters required for a random scan are 

< a, a < 1 and Xf/Xf < 0, d < 1 (74) 

since the diagonal elements of a hermitian matrix must be bounded by its largest and lowest 
eigenvalues. The result of a scan over all these parameters (including the eigenvalues) but 
for fixed A23 is a highly peaked distribution centered at some value. This means that the 
generic value of \Vtb\ is controlled almost exclusively by A23. Of course one may think that 
after fixing the masses to the experimental values, as well as the other mixing angles, maybe 
we might loose this statistical prediction. This is not the case as was shown in Figure [TJ 
where the scan is performed now with only two free parameters, the rest having been fixed 
by other experimental observables. There is still a clear correlation between the value of the 
phase difference and the value of the angle Vtb- 

The imposition of this specific two- angle ansatz (the (13 — 13) case) on the fermion mass 
matrices amounts to 2 constraints per mass matrix, and therefore 4 constraints in the quark 
sector. What we have noted is that with only one more constraint, i.e. the fixing of the 
phase difference A 2 3, we are able to give a statistical prediction for the value of the angle 
I Vtb I (or I Vj-3 1 ) , and this, irrespective of any other parameter. In particular if A 2 3 = we 
would expect \Vtb\ ~ 1 and in the lepton sector, if we fix A 23 — tt/2, then we would expect 
|V T3 | ~ l/y/2, both cases being close to the experimental fits. Two hermitian matrices 
contain 18 free parameters altogether and so it is nontrivial that after imposing only 5 
constraints on the whole set we obtain one prediction, irrespective of anything else. 
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Once noted this nontrivial property, we continued analyzing the rest of consequences of 
our scheme and showed for example how the CP violating phases in both quark and lepton 
sectors depend on the phases of the original mass matrix elements. Another interesting 
outcome was the realization of how to treat in a similar way the quarks and charged lepton 
mass matrices, and use the special case of the neutrino matrix to explain in a transparent 
way the observed differences between the lepton and quark sectors. In particular we also 
analyzed the consequences of further imposing a (2 -h- 3) symmetry in the neutrino sector. 

It is true that out of many possible implementations of the two-angle ansatz we chose to 
study only one case. We felt that this one case would show the main features of the scheme 
in a transparent way, and so we leave the systematic case by case study of the scenario for 
the future, although part of this work is already underway. 
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APPENDIX A: DIAGONALIZATION OF A HERMITIAN MATRIX 

We first note that the diagonal elements of a hermitian matrix are bounded by its smallest 
and largest eigenvalue respectively. This means that in order for a hermitian matrix to have 
zeroes in the diagonal entries, we must have at least one eigenvalue positive and one negative, 
the third one can take either sign, to accommodate the zeroes in the diagonal elements. For 
simplicity, we do not consider this case here and instead concentrate on positive definite 
hermitian matrices such that the eigenvalues are all positive. Of course the results can be 
trivially extended for the case of a more general hermitian matrix, not necessarily definite 
positive. 

Let's introduce our notation by considering the positive definite hermitian matrix H 



( 7 x g\ 



H 



(75) 



x* a b 
\g* b* a) 

The diagonal entries of H must be real and positive and are bounded by its smallest and 
largest eigenvalue respectively. Taking < Ai < A 2 < A 3 the bounds are 

< Ai < 7 < A 3 (76) 
< Ai < a < A 3 (77) 
< Ai < a < A 3 (78) 

The off-diagonal entries of H are also bounded by its smallest and largest eigenvalues in the 
following way; since a 3 x 3 hermitian matrix has only three blocks of off-diagonal elements, 
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namely (x,g), (x,b*) and (g, b) up to conjugation [18]: 



U<VW+W<\\^-M (79) 
0<VW+W<1 |A 3 - Ail (80) 
o< v^T[^< i |A 3 -Ai| (81) 



which means that: 



< |x| < - |A 3 - Ai| (82) 
0<M<^ [A 3 — Aa| (83) 
0< \b\ < l - IA3-A1I (84) 

The above results are valid for any N x N hermitian matrix. These bounds might not be 
very revealing in the quark and charged lepton sectors, since the difference between the 
heaviest and the lightest eigenvalues is of the order of the largest eigenvalue, and so the 
constraint on the off-diagonal entries is quite mild. On the other hand, if the difference 
between the heaviest and lightest eigenvalue is very small then one sees that the off-diagonal 
entries must actually be smaller than half that difference, and so the hermitian matrix must 
be close to diagonal form. This case might be possible in the neutrino sector in the case of 
quasi-degenerate masses. We now write the three invariants Tr(H), Tr(H 2 ) and Det(H): 

Tr(H) = a + a + 7 = Ai + A 2 + A 3 , (85) 
Tr(H 2 ) = 2(\x\ 2 + \b\ 2 + \g\ 2 ) + a 2 + a 2 + 7 2 = A 2 + A 2 + A 2 (86) 
Det(H) = 7(aa - \b\ 2 ) - a\x\ 2 - a\g\ 2 + 2Re(bxg*) = AiA 2 A 3 (87) 

They can be rewritten as 

a = A1 + A2 + A3-C1-7, (88) 
, ,2 = (7 - AQ(A 2 - 7) (A 3 - 7 ) - \g\ 2 (a - 7 ) + 2Re(bxg*) 
X (a - 7) 



(a — Ai)(o — A 2 )(A 3 — a) + \g\ (a — a) — 2Re(bxg 



(90) 



(a -7) 

As noted in the text, the interesting thing of this notation is that the constraint formulae 
on x and b actually become algebraic solutions for both x and b when the term Re(bxg*) 
vanishes identically. 
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When one is interested in finding the eigenvalues of the mass matrix given in Eq. (75), it 
is necessary to find solutions of the characteristic equation: 







(a-A)((A- 7 )(A-a) 



a) — \x\ 



+ (A- 7 )|6| 2 + (X-a)\g\ 2 + 2Re(bxg*) (91) 



In our case, however, we do know the eigenvalues, which are quantities measured experi- 
mentally. Therefore it is preferable to treat them as known parameters instead of unknown 
variables. In this case, we can obtain simple analytical forms for the unitary matrices W u 
and Wd, responsible for diagonalizing the mass matrices H u and Hd- Here is the sim- 
ple procedure: let E u \. = (E^* , E^ , Elf) be one eigenvector of the mass matrix H UJ i.e. 
H u E u x t = Aj E u x t , where Aj is one of the eigenvalues of H u . This means that E,^ = 0, 
or simply 



Det{H uK ) = 



(92) 



which is in fact the characteristic equation of H u given in Eq. (91) and where 

/ (7 - A t ) 



\ 



X 

(a - Xi 
b* 



9 
b 



\ 



(93) 



x 

9* b* (a - Xi) ) 

Because, the Det(H u \ t ) vanishes, we know that one row of the matrix must be a linear 
combination of the other two. Depending on which row we choose to treat as linearly 
dependent, we can obtain different (but equivalent) parametrizations for the eigenvectors 
Eu^. Since the homogeneous equation ( |92| ) can be multiplied by an arbitrary number, we 
will obtain ratios (of 2 x 2 determinants) for the eigenvector components: 
Using the first row as linearly dependent, we obtain: 

x* b 
9* (a-Xi) 



E^/E^ 



E^/Ep 



a — Xi) b 
b* (a - Xi) 

x* (a - Xi) 
9* b* 



(94) 



(a - Xi) b 
b* (a - Xi) 



(95) 
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To simplify the notation, we can choose E t * 
eigenvector 



a — Aj b 
b* a - Xi 



and obtain finally the 



rplst row 



(a - Xi) b 
b* (a - A;) 



x* b 
9* (o-Aj) 



If instead, we choose row 2 as linearly dependent we obtain 



E 



2nd row 



x g 
b* (a - Xi) 



(7 - Aj) g 
9* (a - X, 



x* (a - Xi) 
9* b* 

(7 - Xi) x 
9* b* 



(96) 



(97) 



and finally when the third row is treated as linearly dependent we have 



E 



3d row 



( 


X 


9 




(7 -A 


i) 9 




(7 " \ 


1 X 




(a — X, 


) b 


1 


X* 


b 




X* 


(a — 



(98) 



Since there are three eigenvalues, and three parametrization choices for each, we have 9 
equivalent parametrizations of the unitary matrix W u constructed using any of the previous 
eigenvectors. Note that the difference between each parametrization, is the explicit absence 
of one of the diagonal elements in the formulae (7 in the first one, a in the second one 



and a in the third one). Of course, the characteristic equation (91) relates each 'basis' with 
each other. In the text, we used this parametrization freedom to choose a mixing matrix 
W that has the correct form when b — > 0, g — > and x — > 0. In that limit, the different 
parametrization can be written as: 



771 1st row 
Tp2nd row 
7-i3d row 



((a -Xi)(a-Xi), 0,0) 

(0,-( 7 -A i )(a-A,),0) 

(0,0,(a-A,)( 7 -A l )) 



(99) 
(100) 
(101) 



Some of these, can become problematic since in this limit, we have (7 — Ai) — > 0, (a — A2) — > 
and (a — A3) . There is however, out of the 9 possible combinations, one single 
choice which has the correct smooth asymptotic behavior. The choice is to take row for 



Ai, E. 



W 



2nd row 
u\ 2 



for A 2 and E, 



3d row 



for A 3 , i.e. 



( (a — Xi)(a — Ai) — \b\ 2 gb* — x(a — A2) 

g*b - x*(a - Ai) (7 - A 2 )(a - A 2 ) - \g\ 2 

iVi N 2 N 3 

s*6*-g>-Ai) xg*-b*( 1 -X 2 ) ( 7 -A 3 )(a-A 3 

v m n 2 n 3 



xb - g(a - A 3 ) 
iV 3 

x*g - 6(7 - A 3 ) 



\x\ 



(102) 
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with the normalization parameters 

N 2 = (A 3 - Ai)(A 2 - Ax) [(a - A x )(a - Ai) - |6| 2 ] . 
(A 3 - A 2 )(A 2 - Ax) [(a - A 2 )(A 2 - 7) + \g\ 2 } ■ 
(A 3 - A 2 )(A 3 - Ax) [(A 3 - 7) (A 3 -a)- \x\ 2 ] 



(103) 
(104) 
(105) 



Of course depending on the specific scenario studied, some other parametrization might be 



used. For example, one can use only the eigenvectors of Eq. (96), i.e. 
/ (a- Ai)(o- Ax) - \b\ 2 
g*b — x*(a — Ax) 



W 



x*b* 



iVx 
g*(a 



V 



iVx 



or use the ones from Eq. (97): 
/ gb* - xja - Ax) 



W 



(7 - Ai)(a - Ax) 



\9\ 



\ 



xg 



N'{ 



K 



N'l 



or use the ones from Eq. (98): 
/ xb-g(a-X 1 ) 



W 



x g 



W 
■6(7 



V 



N" 

(7 - Ai)(a - Ax) - \x\ 
Nf 



(a 


- A 2 )(a - A 2 ) - 


\b\ 2 


(a 


- A 3 ) (a - A 3 ) - 


H 2 \ 




AT' 

g*b — x*(a — X 2 ] 






g*b -x*(a- A 3 ) 






AT' 

2 

x*b* — g*{a — A 2 


) 


N J 

x*b* — g*(a — A 3 


) 




N' 
iV 2 






K 


1 

/ 


(7 


gb* — x(a — A 2 ) 
N 2 

- A 2 )(a - A 2 ) - 


\g\ 2 


(7 


gb* - x(a - A 3 ) 

NS 

- A 3 ) (a - A 3 ) - 


\ 

\g\ 2 




N 2 

xg*-b*{ 7 -\ 2 ) 






^3 

xg'-b'h-Xs) 






N 2 






N^ 


J 




xb — g(a — A 2 ) 






xb - g(a - A 3 ) 


\ 




No' 

x*g - 6(7 - A 2 ) 






N 3 

x*g - 6(7 - A 3 ) 




(7 


N? 

- A 2 )(a - A 2 ) - 


\x\ 2 


(7 


N 3 

- A 3 )(a - A 3 ) - 


\x\ 2 










iV 3 


J 



(106) 



(107) 



108) 



where the normalization constants iV", N^', N' 2 , N 2 ', N% and N% can easily be obtained from 
iVx, N 2 and N% by permutation of the mass eigenvalues. 

All in all, one can write the matrix W with 9 different parametrizations (modulo a 



diagonal phase matrix), by permutations of the column vectors from Eqs. (96), (97) and 
(98). Of course, when one takes some special limits some of the parametrizations will reveal 
themselves less useful, as it is possible to find undetermined expressions of the type 0/0. For 



example in the parametrization shown in Eq. ( 107) this will happen if we take simultaneously 
the limits g — > and x — > . In that situation one simply chooses the parametrization 
with a smooth limit. 



It is easy to realize from Eqs. (106), (107) and (108) that one can obtain very simple 
expressions for the absolute value of each element of W. By taking the real elements of these 
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three parametrizations and squaring them we obtain, in terms of the eigenvalues and the 
parameters 7, a, x and g: 

(A 2 -7)(A 3 -7) + M 2 + M 2 



(A 3 -Ai)(A 2 -Ai) 



(109) 



\W 12 \' = 



2 ( 7 -A 1 )(A 3 -7)-k| 2 -|^| 2 



(A 3 -A 2 )(A 2 -A 1 ) 



(110) 



13 



|2 _ N 2 + |^| 2 -( 7 -A 1 )(A 2 - 7 ) 



(A 3 -A 2 )(A 3 -Ai) 



(111) 



, |2 = (7-A 1 )(a-A 1 )-M 2 
1 211 (A 3 -A 1 )(A 2 -A 1 ) 



(112) 



\W 22 r = 



i^| 2 + (A 2 - 7 )(a-A 2 ) 



23| 



(A 3 -A 2 )(A 2 -Ai) 

(A 3 -7)(A 3 -a)-|^| 2 
(A 3 -A 2 )(A 3 -A 1 ) 



(113) 



(114) 
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(7 - Ai)(A 3 + A 2 - 7 - a) - |:r| 2 



32 



\w- 



33 I 





(A 3 -Ax)(A 2 - 


Ai) 


(A 2 


-7)(A 3 + Ai -7 


— a) + x 2 




(A3 — A 2 )(A 2 — 


Ai) 


(A 3 


- 7) (7 + a - Ai - 


-A 2 )-|x| 2 




(A3 — A 2 )(A 3 — 


Ai) 


13) 


ANSATZ 





(115) 



(116) 



(117) 



Let's recall the notation for H u and H d : 







X 


g\ 




(P 


y 


h\ 


H u = 


X* 


a 


b 


, H d = 


y* 


P 


f 




\g* 


b* 


a) 




\h* 


r 


d) 



(118) 



W u and Wd are the unitary transformations diagonalizing H u and Hd respectively. We 
want to find the parametrization of both W u and W d when the elements (W„)i 3 and (W d ) 13 
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vanish. The requirement for the cancellation of these elements is xb — g(a — m t ) — and 
yf — h(j3 — m b ) =0 , and with them we obtain in both sectors the simpler identities: 



\x\ ■ 

\g\ 2 
\b\ 2 



(7 - m u )(m c - 7)(m f - a) 
{2m t — a — a) 

(7 ~ m u )(m c - i){m t - a) 
(2m t — a — a) 

(m t - a)(m t - a) 



1 2/1 



\h\ 



2 _ (p - m d )(m s - p)(m b - (3) 
(2m b -0-d) 

2 = (p - m d ){m s - p)(m b - d) 
(2m b -P-d) 

2 = (m b - P)(m b - d) 



S x -5 g + 5 b 



IT 



Sy - 5 h + 5 



f 



IT 



(119) 



From the above expressions we obtain the following simple parametrization of up and down 
quark mixing matrices, with (W u )i3 = 0, (Wd)i3 = : 



(m c -7) 
(m c -m u ) 



il-m u ) ( 
(m c -m u ) 





(■y-m u )(m t -a) ^i5 x 
(m c -m u )(2m t -a-a) 

(m c -j)(mt-a) 
(m c -m u )(2m t -a-a) 

(mt-a) -jS b 



(n-m u ){m t -a) c i(5 x +8 h ) \ 
(m c -j)(m t -a) „i5 h 



(m c -m u ) (2m t -a-a) 



(2mt—a—a) 



(mt—a) 
(2mt—a—a) 



(120) 



and 



W,= 



( 



\ 



(m s -p) 
(m s -m d ) 



(p-™d) c i$ v 
(m s -rn d ) 



(p-m d )(m b - /3) 



(m s -m d )(2m b -l3-d) 



(p-m d )(m b -d) 



-i(5y+S f ) 



(m s -p)(m b -0) 
{m a -m d )(2m b - f3-d) 

e" 



e io f 



(m s -p)(m b -d) p-i5f 
(m s -m, d )(2m b —f)—d) 



(m b -d) 
(2m b -P-d) 

(m b -/3) 
(2m b -/3-d) 



(121) 



/ 



(m s -m d ){2m b -P-d) 

One can now check that in this ansatz, which basically shows that the structure of the quark 
mixing matrix is 

Ci-y S p 6 ^ I S fy C p 6 S sy S 63 N 



/ c 7 c p + SpS^e^-^C* 

SsyCp (3 | C-y S pG ^ (^J 

V 



-s p S*e~ i6 y 



-c p S* 



(122) 



C J 



where C, S, c 7 and c p are given by 

(y(m t -a)(m b -d) e^f'^ + ^(m* - a)(m b - /?)) 
^/ (2m t — a — a) (2m b — (5 — d) 
(m t — a)(m b — (3) — e l ^ 5 f~ 5b ^yj{m t — a)(m b — d)j 



C = 



V, 



tb 



s = 



■sj (2m t — a — a) (2m b — (5 — d) 



1 m c 


-7 


m c - 


- m u 


1 m s 


- P 



m s - m d 



(123) 

(124) 

(125) 
(126) 
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and Si = a/1 — cf and |C| 2 + |S'| 2 = 1. The previous form for Vckm can then be put in the 



form given in Eq.(28) by pulling out two diagonal phase matrices P u and Pd as given in the 
main text. 

We can also quickly compute the Jarlskog invariant J for this mixing matrix, for example 
by computing 

J = im(v cb v t tv; d v td ) 

J = Im (-c*s 2 p \SC\ 2 + GyCpSySplSfe^-^C*) 

( y/(m t -a)(m b -d) Sm (S b ~5 j+5 x S y ) + y/(m t -a)(m b -P) Sm (S x S y ) ) 

J = c lWp \S\ 2 ± = L (127) 

\J {2mt—OL—a)(2m b —fi—(£) 
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